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A classical view of neural coding relies on temporal firing synchrony among functional 
groups of neurons, however, the underlying mechanism remains an enigma. Here we 
experimentally demonstrate a mechanism where time-lags among neuronal spiking leap 
from several tens of milliseconds to nearly zero-lag synchrony. It also allows sudden 
leaps out of synchrony, hence forming short epochs of synchrony. Our results are 
based on an experimental procedure where conditioned stimulations were enforced 
on circuits of neurons embedded within a large-scale network of cortical cells in vitro 
and are corroborated by simulations of neuronal populations. The underlying biological 
mechanisms are the unavoidable increase of the neuronal response latency to ongoing 
stimulations and temporal or spatial summation required to generate evoked spikes. These 
sudden leaps in and out of synchrony may be accompanied by multiplications of the 
neuronal firing frequency, hence offering reliable information-bearing indicators which may 
bridge between the two principal neuronal coding paradigms. 
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INTRODUCTION 

One of the major challenges of modern neuroscience is to eluci- 
date the brain mechanisms that underlie firing synchrony among 
neurons. Such spike correlations with differing degrees of tempo- 
ral precision have been observed in various sensory cortical areas, 
in particular in the visual (Eckhorn et al., 1988; Gray et al., 1989), 
auditory (Ahissar et al., 1992; Nicolelis et al., 1995), somatosen- 
sory (Nicolelis et al., 1995), and frontal (Vaadia et al., 1995) areas. 
Several mechanisms have been suggested, including the slow and 
limited increase in neuronal response latency per evoked spike 
(Vardi et al., 2013b). On a neuronal circuit level its accumulative 
effect serves as a non-uniform gradual stretching of the effective 
neuronal circuit delay loops. Consequently, small mismatches of 
only a few milliseconds among firing times of neurons can vanish 
in a very slow gradual process consisting of hundreds of evoked 
spikes per neuron. 

The phenomenon of sudden leaps from firing mismatches of 
several tens of milliseconds to nearly zero-lag synchronization, 
below a millisecond, is counterintuitive. Since the dynamical vari- 
ations in neuronal features, e.g., the increase in neuronal response 
latencies per evoked spike, are extremely small, one might expect 
only very slow variations in firing timings. Moreover, relative 
changes among firing times of neurons require dynamic relax- 
ation of the entire neuronal circuit to achieve synchronization. 
Hence, sudden leaps, in and out of synchrony, seem unexpected. 

In the present study, we propose a new experimentally cor- 
roborated mechanism allowing leaps in and out of synchrony. 
The procedure is based on conditioned stimulations enforced 
on neuronal circuits embedded within a large-scale network of 



cortical cells in vitro (Marom and Shahaf, 2002; Morin et al, 
2005; Wagenaar et al., 2006; Vardi et al., 2012). These stimulations 
varied in strength, so that the evoked spikes of selected neurons 
required temporal summation. We demonstrate that the underly- 
ing biological mechanism to sudden leaps in and out of synchrony 
is the unavoidable increase of the neuronal response latency 
(Aston-Jones et al, 1980; De Col et al, 2008; Ballo and Bucher, 
2009; Gal et al., 2010) to ongoing stimulations, which imposes a 
non-uniform stretching of the neuronal circuit delay loops. 

MATERIALS AND METHODS 
CULTURE PREPARATION 

Cortical neurons were obtained from newborn rats (Sprague- 
Dawley) within 48 h after birth using mechanical and enzymatic 
procedures (Marom and Shahaf, 2002; Vardi et al, 2012, 2013b). 
All procedures were in accordance with the National Institutes of 
Health Guide for the Care and Use of Laboratory Animals and 
Bar-llan University Guidelines for the Use and Care of Laboratory 
Animals in Research and were approved and supervised by the 
Institutional Animal Care and Use Committee. 

The cortex tissue was digested enzymatically with 0.05% 
trypsin solution in phosphate-buffered saline (Dulbecco's PBS) 
free of calcium and magnesium, supplemented with 20 mM 
glucose, at 37°C. Enzyme treatment was terminated using 
heat-inactivated horse serum, and cells were then mechanically 
dissociated. The neurons were plated directly onto substrate- 
integrated multi-electrode arrays (MEAs) and allowed to develop 
functionally and structurally mature networks over a time period 
of 2-3 weeks in vitro, prior to the experiments. Variability in 



Frontiers in Neural Circuits www.frontiersin.org October 2013 | Volume 7 | Article 176 | 1 



Vardi et al. 



Sudden synchrony leaps 



the number of cultured days in this range had no effect on 
the observed resuhs. The number of plated neurons in a typ- 
ical network is in the order of 1,300,000, covering an area of 
about 380 mm^. The preparations were bathed in minimal essen- 
tial medium (MEM-Earle, Earle's Salt Base without L-Glutamine) 
supplemented with heat-inactivated horse serum (5%), glu- 
tamine (0.5 mM), glucose (20 mM), and gentamicin (lOg/ml), 
and maintained in an atmosphere of 37° C, 5% CO2 and 95% air 
in an incubator as well as during the electrophysiological mea- 
surements. All experiments were conducted on cultured cortical 
neurons that were functionally isolated from their network by a 
pharmacological block of glutamatergic and GABAergic synapses. 
For each plate, 12-20 [lI of a cocktail of synaptic blockers was 
used, consisting of 10 (jlM CNQX (6-cyano-7-nitroquinoxaline- 
2,3-dione), 80 |iM APV (amino-5-phosphonovaleric acid), and 
5 |xM BicucuUine. This cocktail did not block the spontaneous 
network activity completely, but rather made it sparse. At least 
1 h was allowed for stabilization of the effect. 

MEASUREMENTS AND STIMULATION 

An array of 60 Ti/Au/TiN extracellular electrodes, 30 [im in diam- 
eter and spaced either 200 or 500 |xm from each other (Multi- 
Channel Systems, Reutlingen, Germany) was used. The insulation 
layer (silicon nitride) was pre-treated with polyethyleneimine 
(Sigma, 0.01% in 0.1 M Borate buffer solution). A commer- 
cial setup (MEA2100-2x60-headstage, MEA2100-interface board, 
MCS, Reutlingen, Germany) for recording and analyzing data 
from two 60-electrode MEAs was used, with integrated data 
acquisition from 120 MEA electrodes and 8 additional analog 
channels, integrated filter amplifier and 6-channel current or volt- 
age stimulus generator (for both MEAs). Mono-phasic square 
voltage pulses (—900 to — lOOmV, 100-500 (is) were applied 
through extracellular electrodes. Each channel was sampled at a 
frequency of 50 k sample/s. Action potentials were detected on- 
line by threshold crossing. For each of the recording channels a 
threshold for spike detection was defined separately, prior to the 
beginning of the experiment. 

CELL SELECTION 

Each circuit node was represented by a stimulation source (source 
electrode) and a target for the stimulation — the recording elec- 
trode (target electrode). These electrodes (source and target) were 
selected as the ones that evoked well-isolated, well-formed spikes 
and reliable responses with high signal-to-noise ratio. This exam- 
ination was done with stimulus intensity of — SOOmV using 30 
repetitions at a rate of 5 Hz followed by 1200 repetitions at a rate 
of 10 Hz. 

STIMULATION CONTROL 

A node response was defined as a spike occurring within a typi- 
cal time window of 2-10 ms following the electrical stimulation. 
The activity of all source and target electrodes was collected, 
and entailed stimuli were delivered in accordance to the circuit 
connectivity. 

Circuit connectivity, t 

Conditioned stimulations were enforced on the circuit neurons 
embedded within a large-scale network of cortical cells in vitro. 



according to the circuit connectivity. Initially, each delay was 
defined as the expected time between the evoked spikes of two 
linked neurons; e.g., conditioned to a spike recorded in the target 
electrode assigned to neuron A, a spike will be detected in the tar- 
get electrode of neuron B after xab ms. For this end, conditioned 
to a spike recorded in the target electrode of neuron A, a stimu- 
lus will be applied after xab-Lb{0) ms to the source electrode of 
neuron B, where LsiO) is the initial latency of neuron B. 

In cases where missed evoked spikes caused a termination of 
the neuronal circuit activity, stimulation was given to neuron A 
after a period of 100 ms, to restart the circuit's activity. 

All neurons were stimulated at a rate of 10 Hz (Figures 1, 3) or 
8 Hz (Figure 2), before the leap to synchronization. 

Strong stimulations, (— SOOmV, 200 (jls), resulting in a reli- 
able neural response, were given to all circuit neurons exclud- 
ing neuron C (Figures 1, 2) and E (Figures). Weak stimu- 
lations (Figure 1: — 450mV, 40 (xs. Figure 2: — 600mV, 60 jis. 
Figure 3: — 700mV, 60(jls) were given to neuron C (Figures 1, 
2) or E (Figure 3), so that an evoked spike is expected 
only if the time-lag between two consecutive weak stimula- 
tions is short enough. In cases where the time-lag between 
two consecutive stimulations was shorter than 20 [is (from 
the end of the first stimulation to the beginning of the con- 
secutive one), a unified strong stimulation was applied, to 
overcome technical limitations. The weak stimulations were 
defined for each neuron separately, due to differences in their 
threshold. 

Tjs (TS stands for temporal summation) is the maximal time- 
lag between two weak stimulations which typically results in an 
evoked spike. This quantity was empirically estimated by gradu- 
ally changing the time-lag between two weak stimulations, and 
found to differ between neurons. 

DATA ANALYSIS 

Analyses were performed in a Matlab environment (Math Works, 
Natwick, MA, USA). Action potentials were detected by threshold 
crossing. In the context of this study, no significant difference was 
observed in the results under threshold crossing or voltage min- 
ima for spike detection. Reported results were confirmed based 
on at least ten experiments each, using different sets of neurons 
and several tissue cultures. 

RESULTS 

LEAP TO SYNCHRONY ACCOMPANIED BY A DOUBLED FIRING 

FREQUENCY 

Experimental results 

We first demonstrate leaps to synchrony using a neuronal cir- 
cuit consisting of four neurons and conditioned stimulations 
split into weak/strong stimulations (Figure lA). A strong stim- 
ulation consists of a relatively high amplitude and/or relatively 
long pulse duration such that an evoked spike is generated reli- 
ably, whereas a weak stimulation consists of a lower amplitude 
and/or pulse duration, such that an evoked spike is expected only 
if the time-lag between two consecutive weak stimulations is short 
enough. All delays (denoted on connecting lines between neu- 
rons in Figure lA) were selected to initially include the response 
latency of the target neuron, e.g., the time-lag from neuron A 
to B, Xab, was initially set to x-Lb{0) where -Lb(O) stands for 
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FIGURE 1 I A sudden leap to synchrony accompanied by frequency 
doubling. Notations used: Sync^g, the absolute time-lag between the spikes 
of neurons A and B; A(Stimc), the absolute time difference between two 
weak stimulations to neuron C; ALq, the increase in response latency of 
neuron D after n evoked spikes. (A) Schematic of a neuronal circuit consisting 
of four neurons and weak/strong stimulations represented by dashed 
(green)/full (black) lines. An initial stimulation is given to neuron A. (B) 
Experimental measurements of A(Stimc) as a function of the spikes of 
neuron A. AlStim^) is initially set to e aj 0.8 ms (green line) with x = 50 ms 
and Tfs 0.24 ms (presented by the dashed horizontal green line). A unified 
longer stimulation was given in events where the time-lag between the weak 
stimulations <20|ls [presented by A(Stimc) = 0]. Sync^g is presented by the 
blue line, indicating a sudden leap from x = 50 ms to nearly zero-lag 
synchronization. (C) Spike trains of the four neurons. A sudden leap to 
Sync/\B 0 occurs at time/2x = 122.5 (at spike 121 of neuron A) immediately 
following a single evoked spike of neuron C. It is accompanied by a doubled 



firing frequency, from ~10 to ~20 Hz. Sync^g 0 is robust to response 
failures of neuron C, e.g., time/2x = 124.5. (D) Sync^g as a function of the 
spikes of neuron A, for various e, where the data for e = 0.8 (blue) is the 
same as in (B,C). The number of spikes to a leap to synchrony increases with 
e. (E) ALo for repeated stimulations at 10 Hz. ALo at the synchrony leap for 
different e are colored following (D). Note that Spikeo is equal to Spike^ in 
(B,D). (F) Results of population dynamic simulations where each neuron in 
(A) is now represented by a population comprised of 40 Hodgkin-Huxley 
neurons, each one innervated by four randomly chosen neurons from each of 
its driving clusters. The delays between neurons are taken from a Gaussian 
distribution centered at the delays of the single neuron case with a variance 
of 0.2 ms. For simplicity, each time a neuron fires all of its outgoing delays are 
increased by 0.04 ms. The simulation parameters were e = 2 ms and 
Tfs 1 .3 ms. (G) Raster plot of the 120 neurons comprising nodes A, B, and 
C. A leap to synchrony occurs at time/2x 20, accompanied by a doubling of 
the firing frequency. 



the initial response latency of neuron B. For x = 50 ms, neu- 
rons A and B initially fire alternately, in and out of phase, at 
a frequency of ~10Hz (Figure IB). Neuron D fires ~t/2ms 
laggard to neuron A (Figure IC) and the time-gap between 



two weak stimulations arriving at neuron C, A(Stimc), is ini- 
tially g (Figures 1A,B). The experimentally estimated maximal 
time-gap between stimulations of neuron C which generates an 
evoked spike (temporal summation) is denoted by Tj-g, thus for 
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FIGURE 2 I A sudden leap to synchrony accompanied by tripled 
frequency. Notations used: Sync/ie, tlie absolute time-lag between the 
spikes of neurons A and B; A(Stimc), the absolute time difference 
between two weak stimulations to neuron C; ALq, the increase in 
response latency of neuron D after n evoked spikes. (A) Schematic of a 
neuronal circuit as in Figure 1A, however, the delay from neuron B to A 
is now 2x. (B) Experimental measurements of A(Stimc), similar to 
Figure IB, with B«i 0.5ms, 3x = 125 ms and Tjs^O.l (presented by the 
dashed horizontal green line). Sync^g. (blue line) indicating a sudden leap 
from T as 125/3ms to nearly zero-lag synchronization. (C) Spike trains of 
the four neurons. A sudden leap to synchronization, Sync,4B ^ 0, occurs 
at time/3x = 44 (at spike 44 of neuron A) consecutive to three evoked 
spikes of neuron C. This is accompanied by tripled firing frequency of 
neurons A and B, from ~8 to ~24 Hz. Sync^ig 0 is robust to response 



failures of neuron C, e.g., time/3x = 46.33. (D) Sync^g 3S a function of 
the spikes of neuron A for various e, where the number of spikes to the 
leap to synchrony increases with e. The data for s=0.5 (blue) is the 
same as in (B,C). The observed oscillations in Sync^g before a leap to 
synchrony originate from response failures of neuron C, and similarly 
oscillations in a leap out of synchrony originate from response failure of 
either neuron A or B. (E) ALq, for repeated stimulations at 8 Hz. ALg at 
the leap for different e are indicated and colored following (D), 
approximately verifying Equation (1), e.g., for e = 0.8ms and Tjs 0.2ms, 
A/.d(197) gives ~0.6ms. Note that Spikeo is equal to Spike^ in (B,D). 
(F) Results of population dynamic simulations similar to Figures 1F,G 
with E = 2ms, 7rs!« 1.3ms and 3x= 125ms. (G) Raster plot of the 120 
neurons comprising nodes A, B, and C. A leap to synchrony occurs at 
time/3x «i 20, accompanied by tripled firing frequency. 



A(Stimc) > Tjs ^ 0.24 ms neuron C typically does not fire. As 
a result of the increase in the response latency of neuron D, 
A(Stimc) is reduced (green line in Figure IB) sufficiently so that 
neuron C starts firing [A(Stimc) < Tjs] (Figure IC). The circuit 
now consists of two delay loops, ~2t (A-B-A) and ~3t (A-C-B- 
A). Since the greatest common divisor (GCD) of the circuit delay 
loops is GCD(2,3) = 1, conditioned to the firing of neuron C, 



zero-lag synchronization between neurons A and B is theoreti- 
cally expected (Kanter et al., 2011) after a very short transient, x 
(Figure IC). This phenomenon is quantitatively measured by the 
time-lag between spikes of neurons A and B, Sync^^. The emer- 
gence of zero-lag synchrony is clearly demonstrated by the leap 
to SyncAB ~ 0 ms (blue line in Figure IB), which is accompa- 
nied by a sudden frequency multiplication from ~10 to ~20Hz 



Frontiers in Neural Circuits 



www.frontiersin.org 



October 2013 | Volume 7 [ Article 176 | 4 



Vardi et al. 



Sudden synchrony leaps 




FIGURE 3 I Short epochs of synchrony not accompanied by a change 
in frequency. Notations used: SyncAE. the absolute time-lag between the 
spikes of neurons A and E; A(Stimf), the absolute time difference 
between two weak stimulations to neuron E; AL, defined as 
ALb + ALc + ALo — ALf. (A) Schematic of a neuronal circuit consisting of 
six neurons and weak/strong stimulations represented by dashed 
(green)/full (black) lines. (B) Experimental measurements of AfStimg), 
similar to AfStimc) in Figure IB, with e 1.7 ms, x = 50ms and 
775 ^ 0.5 ms (presented by the dashed horizontal green line). The time 
delay between neurons A and E, ~2t, is denoted by the dashed horizontal 
black line. The firing region of neuron E (blue dots bounded by dashed 
vertical guidelines), which is at nearly zero-lag synchronization with the 
firing of neuron A, Sync/i£ »i 0, starts after 77 spikes of neuron A. The 



temporary firing of E terminates after ~200 spikes of neuron A. (C) Spike 
trains of neurons A, E and E, indicating a steady firing frequency (~10Hz) 
of the neuronal circuit independent of the firing of neuron E, where an 
epoch of synchrony, Sync/if 0, begins at time/2T = 77 (at spike 77 of 
neuron A). (D) The number of spikes prior to the firing of neuron E 
increases with e. The mild increase in the firing mismatch, Sync^e, is 
attributed to the additional increase by s of the initial 2x delay loop (E fires 
~2x-|-E laggard to A, however, the time-gap between consecutive firings 
of A is ~2x-|-2e). The data for e= 1.7 (blue) is the same as in (B,C). (E) 
AL for repeated stimulations at 10 Hz. AL at the synchrony leap for 
different e are colored following (D). The number of spikes per neuron 
(e.g., Spike^), n, until the leap to synchrony increases with b and can be 
obtained from Equation (1), where ALq is substituted by AL. 



(Figure IC). Note that the sudden multipUcation in frequency, by 
itself, shortens SyncAB from 100 to 50 ms, however, it cannot lead 
to zero-lag synchrony. The sudden emergence of Sync^B ~ 0 ms 
requires only a single firing of neuron C, and is then main- 
tained by the mutual firing of neurons A and B, independently 
of the firing of neuron C (Figure IC). For a given Tjs, the num- 
ber of evoked spikes of neuron D until the leap to synchrony, 
«, increases with e (Figure ID). Quantitatively, using the exper- 
imental response latency profile of neuron D, Ld, one can find n 
fulfilling the equality: 

ALD(n)^e-TTs (D 

where ALd(m) stands for the increase in response latency of neu- 
ron D after n evoked spikes (Figure IE). Note that neuron D is 
laggard to neuron A, thus the number of evoked spikes of neuron 
A until the leap to synchrony increases with e as well, in accor- 
dance with Equation (1) (Figure ID). Since Tjs varies between 
neurons and even within the same neuron over different trials. 



deviations from this equation are expected (e.g., ALd for s = 
0.8 ms and 6 = 1 ms are almost the same. Figures 1D,E). A slow 
gradual increase in Sync^B after a leap to synchrony (Figure ID) 
is theoretically attributed to the difference in the increase of 
neuronal response latencies |ALa(m)-ALb(«)| and the leap out 
of synchrony (Figure ID) is a consequence of a response fail- 
ure of neurons A and/or B (see Section "Slow Divergence out 
of Synchrony" in Appendix). Similar results were obtained and 
exemplified for spatial summation (not shown), where weak stim- 
ulations were given to a neuron through two different source elec- 
trodes. An evoked spike is expected only if the time-lag between 
two consecutive weak stimulations, controlled by the relative 
stimulation timings of the source electrodes, is short enough. 

Simulations of population dynamics 

The sudden leap to synchrony was experimentally verified under 
the limitation where each circuit node is represented by a sin- 
gle neuron, and is demonstrated to be robust under simulations 
of population dynamics (Figures 1F,G). Each one of the four 
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nodes (Figure lA) now represents a population comprised of 
40 Hodgkin-Huxley sparsely connected neurons (for simulation 
details, see Vardi et al., 2013a). For the parameters used, Tjs ~ 
1.3ms, 8 = 2ms and 0.2ms variance for the Gaussian distri- 
bution of the delays, a leap to synchrony is expected following 
Equation (1) after ~20 spikes of cluster A (Figure IF). The sim- 
ulated SjncAB is defined as the absolute difference between the 
average spiking times of the neurons comprising clusters A and B, 
where at least 50% of the neurons in a cluster fired (Figure IG). 
Initially, several neurons in cluster C fire as a result of relatively 
close stimulations from either cluster A or D. This sporadic fir- 
ing is a consequence of the Gaussian distribution of the delays 
between populations, however, their impact on the firing activity 
of cluster B is negligible. As neurons of cluster D fire repeat- 
edly, A(Stimc) decreases and more neurons from cluster C fire. 
Consequently, the activity of cluster C is enhanced such that a leap 
to synchrony is observed, accompanied by frequency doubling 
from ~10 to ~20 Hz (Figures 1F,G). A leap out of synchrony was 
not observed in the simulations, since population dynamics are 
more robust to a single neuron's response failure in comparison 
to a neuronal circuit where each node is represented by a sin- 
gle neuron (Figures 1A,D). Low connectivity, as well as a wider 
Gaussian distribution of delays between populations are expected 
to enhance fluctuations and response failures, and will eventually 
lead to a leap out of synchrony. 

Population dynamics exhibit consistency with most of the 
experimental results, hence minimizing the possibility of these 
results as being only an artifact of the tissue culture. Nevertheless, 
the verification of our results in more realistic scenarios is 
required, including shorter delays and their interplay with the 
neuronal refractory period, the morphology of the neurons 
instead of considering neurons as points (Doiron et al., 2006), 
as well as possible adaptation mechanisms in the form of short 
and long term synaptic plasticity (Abbott and Regehr, 2004; 
Izhikevich, 2006). In addition, more accurate and systematic sta- 
tistical measures of synchrony (Kreuz et al., 2009; Shimokawa and 
Shinomoto, 2009) can be adopted to describe the transition to 
synchrony in the case of population dynamics. 

LEAP TO SYNCHRONY ACCOMPANIED BY TRIPLED FIRING FREQUENCY 

More general features of a sudden leap to synchrony are exem- 
plified by increasing the delay from neuron B to A, tbAj from x 
(Figure lA) to 2x (Figure 2A). The circuit now consists of two 
delay loops, ~3t (A-B-A) and~4T (A-C-B-A) (Figure 2A). Since 
GCD(4,3) = 1, zero-lag synchronization is theoretically expected, 
conditioned to the firing of neuron C. Initially, neurons A and 
B fire at a frequency of ~8 Hz (3t = 125 ms) (Figure 2C) and 
SjncAB ~ 1 (Figure 2B). Neuron C starts to fire as A(Stimc) < 
Tts ^ 0.2 ms, resulting in Sync^iB ^ 0 which is accompanied 
by tripled firing frequency (Figure 2C). The number of evoked 
spikes by neuron D (or its leader neuron A) to the leap increases 
with E in a non-linear manner following ALo{n), in accordance 
with Equation (1) (Figures 2D,E). 

Typically, several leaps in and out of synchrony between neu- 
rons A and B occur before arriving at a stable nearly zero-lag 
synchronization (Figure 2D). These oscillations are attributed to 
unreliable responses of neuron C, and increase the duration of 



the relaxation to synchrony (Figure 2D). Similar oscillations on 
the way out of synchrony (Figure 2D) are attributed to the first 
response failure of either neuron A or B. Consequently, neurons 
A and B fire alternately in time-lags t and 2x. The final exit out 
of synchrony occurs in the second response failure of neurons 
A or B. 

Simulation results (Figures 2F,G) confirmed the robustness 
of the experimentally observed leap to synchrony in popula- 
tion dynamics. The oscillations in the relaxation to synchrony 
are attributed to response failures of cluster C. These failures 
are a consequence of fluctuations in the firing timings of clus- 
ters A and D and the Gaussian distribution of their delays to 
cluster C. 

EPOCHS OF SYNCHRONY NOT ACCOMPANIED BY A CHANGE IN 
FREQUENCY 

A mechanism to leap out of synchrony as well as the inter- 
relation between the sudden leap to synchrony and the firing 
frequency are at the center of the next examined neuronal circuit 
(Figure 3A). This circuit consists solely of a 2t-delay loop, hence 
neurons A and F fire alternately in ~t ms time-lags. Nevertheless, 
neuron A affects neuron E by weak stimulations arriving from 
two comparable initial delay routes; ~2t ms (A-F-E) and ~2t- 
ems (A-B-C-D-E) (Figure 3A). Initially, neuron E does not fire 
since e ~ 1.7 ms > Tjs ~ 0.5 ms. Since the overall increase in 
the neuronal response latency of a chain is accumulative, propor- 
tional to the number of neurons it comprises, A (Stim^) gradually 
decreases below Tjs (Figure 3B) and neuron E suddenly starts to 
fire. Consequently, since neuron A fires every ~2t ms and neu- 
ron E fires ~2t ms laggard to A, Sync^E ~ 0 (Figures 3B,C). As 
A (Stim^) decreases, the response of neuron E becomes more reli- 
able (Figures 3B,C) and a leap out of synchrony is observed when 
A(Stim£) again exceeds ^Tjs (Figure 3B). Since neuron E's fir- 
ing does not close a new neuronal loop, the leaps in and out of 
synchrony do not affect the firing frequency of the neuronal cir- 
cuit (Figure 3C). The number of spikes to synchrony increases 
with E as well as the time-gap between neurons during synchro- 
nization, SyncAB (Figures 3D,E)- Simulation results (not shown) 
confirmed the robustness of the experimentally observed leap in 
and out of synchrony without a frequency change in population 
dynamics. 

DISCUSSION 

Understanding the brain mechanisms that underlie firing syn- 
chrony is one of the great challenges of neuroscience. There are 
many variants of population codes, where a set of neurons in 
a population acts together to perform a specific computational 
task (Palm, 1990; Eichenbaum, 1993; Ainsworth et al, 2012). 
There is much discussion over whether rate coding or temporal 
coding is used to represent perceptual entities in populations of 
neurons in the cortex. A number of reports suggest that almost 
all the information in a stimulus is embedded in the rate code 
of active neurons (Aggelopoulos et al., 2005), while others sug- 
gest that synchrony among spiking of neuronal populations carry 
the information (deCharms and Merzenich, 1996). Experimental 
support for changes solely in firing rate when the perceptual 
task is modified (e.g., Lamme and Spekreijse, 1998; Roelfsema 
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et al, 2004) is as compelling as those works that show changes in 
synchrony in the absence of firing rate changes (e.g., Womelsdorf 
et al, 2005), whereas in other experiments changes in both rate 
and spike correlations are observed concurrently (e.g., Biederlack 
et al, 2006). In any case, the usefulness of rate coding and tem- 
poral coding as information carriers of brain activity is a function 
of the decoding complexity, which is tightly correlated with their 
accuracy. 

Rate and temporal coding are typically inaccurate in brain 
activities, although there are several well-known exceptions where 
neurons fire with high temporal accuracy (Bullock, 1970; Bullock 
et al., 1972; Moortgat et al., 2000). Rate precision, measured by 
inter-spike interval (ISI) distributions, typically follows a broad 
distribution, deviating from a Poissonian one (Amarasingham 
et al., 2006). In addition, relative spike timings between coactive 
neurons are usually within the precision of several milliseconds 
(Kayser et al., 2010; Wang, 2010). In the case of a broad dis- 
tribution of ISIs, the mission to grasp gradual changes in tem- 
poral and/or rate coding (e.g., changes from an average firing 
rate of 5-6 Hz), on a timescale of a few ISIs, is a heavy com- 
putational mission which might not be satisfactorily resolved. 
The underlying cause of this computational difficulty is the 
broad distribution of the ISIs which is overlapped between 
gradually changed temporal codes or gradually changed rate 
codes. 

To overcome this difficulty we proposed a mechanism which 
enables the emergence of a sudden leap to synchrony together 
with or independent of a leap in the firing frequency. This mech- 
anism results in leaps from firing mismatches of several dozens 
of milliseconds to nearly zero-lag synchronization, and can be 
accompanied by a sudden frequency multiplication of the neu- 
ronal firing rate. These sudden changes occur on a time scale 
of extremely few ISIs, and are easUy detectable as the distribu- 
tions of the ISIs before and after the leaps are non-overlapping. 
Hence, one ISI is sufficient to detect the transition without 
accumulatively estimating the ISI distribution. These fast and 
robust indicators might be used as reliable information carriers 
of time-dependent brain activity. 

The proposed mechanism also allows for the simultaneous 
emergence of sudden leaps in rate and temporal synchrony, hence 
bridging between these two major schools of thought in neuro- 
science (Eckhorn et al., 1988; Gray et al, 1989; Ahissar et al., 
1992; Nicolelis et al., 1995). This mechanism requires recurrent 
neuronal circuits, and synchrony appears even among neurons 
which do not share a common drive. Sub-threshold stimula- 
tions (e.g., the stimulations to neuron C in Figures 1, 2 and 
to neuron E in Figure 3) serve as a switch that momentarily 
closes or opens a loop in the neuronal circuit. The state of the 
switch changes a global quantity of the network, the GCD of 
the entire circuit's loops, which determines the state of syn- 
chrony (e.g., zero-lag synchrony, cluster synchrony, shifted zero- 
lag synchrony) (Kanter et al., 2011; Nixon et al., 2012). These 
demonstrated prototypical examples call for a theoretical exam- 
ination of more structured scenarios, including multiple leaps 
in and out of synchrony. In addition, a more realistic biological 
environment has to be examined containing synaptic noise and 
adaptation. 
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APPENDIX 

SLOW DIVERGENCE OUT OF SYNCHRONY 

The slow increase in SynCyiB (Figure ID) is analytically exam- 
ined below for a case of two phase-to-phase neurons, A and B, 
as depicted in Figure Al. The derivation below is in the spirit of 
Ermentrout's analysis of coupled type I membranes (Ermentrout, 
1996). We first define the following quantities and assumptions: 

tiiq) = the timing of the q*^ spike of neuron i, e.g., fA(0) is the 
timing of the first spike of neuron A, where the count starts at 0. 

Liiq) = neuronal latency of neuron i at its spike. 

The initial time delays are tab = '^ba = t- 

Assuming initial conditions, f = 0, where both neurons fire 
simultaneously, i.e., ^4(0) = 0, tsiO) = 0. 

The spiking times of neurons A and B are given by 



(i) tsiq) = tA(q-l) + X + Leiq) 
(ii) tA(q) = tsiq - 1) + t + LAiq) 




Spike 



Spike 



150 200 



FIGURE Al I Slow divergence out of synchrony between two 
phase-to-phase neurons. Notation used: Sync/tg, the time-lag between 
the spikes of neurons A and B. (A) Schematic of two bidirectional 
interconnected spiking neurons. The initial delays between the neurons are 
equal, tab = te/i = t- (B) Response latency of both neurons as a function 
of spike number. The latencies were taken to be = 0.5 x ln[q+ 5) + 2, 
Lb = 0.3 X sqrt((7 -I- 2) + 3, qualitatively similar to latency profiles observed 
in experiments. (C) Sync^g a function of spike number for the latencies 
depicted in (B), assuming Sync/iB(0) = 0. The calculation (brown line) was 
done using Equation (Al) and is in a good agreement with straightforward 
simulations of exact spike times (black dots). For simplicity, the simulated 
SynCy^B is only displayed for even numbers of spikes. 



Substituting (ii) into (i) and vice versa: 

I tsiq) = tsiq - 2) + T + LAiq - 1) + t + Lsiq) 
\ tA(q) = tA(q - 2) + t + Leiq - 1) + x + LAiq) 

one can find that the solution of these coupled recursive equations 
is given by: 



tsiq) = Ej^i Lb {2q') + La {2q' - l) + 2t 
tAi^) = E| = 1 La {2q') + Lb {2q' - l) + 2t 



Consequently, the firing time-gap between the two neurons is 
given by 

SyncAB(<?) = \tBiq) - tAiq)\ 
SyncABiq) = Yll^i + La (2<j' - l) 

- La {2q') - Lb {2q' - l) 
SyncABiq) = J^l^^ {La {2q' - l) - La {2q')) 

-E} = i(^B(2q'-l)-LB(2g')) 

Under the assumption of continuous increase in latency and 
large q 



dL, {2q') 
d(2q') 

SyncABiq) 



(I, {24 - 1) - L, {2q')) > 0 



i 



dLA(2^^ 
d(2q') d(2q') ^ ^' 



rtdM2^ dq' 
Jo d(2q') d(2q') ^ ^' 

|0.5 {LAiq)-LAiQ))-Q.5 {LBiq)-LB (0))| (Al) 



SjncABiq) 



Note that these calculations refer to even values of q. Similar equa- 
tions can be obtained for odd values of q (not shown) . In addition, 
fluctuations in the latencies may also enhance the deviation from 
synchronization. 
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